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Abstract 

We perform a QCD sum rule calculation to determine the mass and the parity of the 
lowest lying pentaquark state. We include operators up to dimension d = 13 in the OPE 
and the direct instanton contributions. We find evidence for a positive parity state. The 
contribution from operators of dimension d > 5 is instrumental in determining the parity 
of the state and achieving the convergence of the sum rule. 
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1 Introduction 



The experimental and theoretical status of 0^-pentaquark remains controversial PP, j2]. 
The QCD sum rules (SRs) have shown to be a very powerful tool for the investigation 
of the properties of conventional and exotic multiquark hadronic states Several 
attempts to describe the properties of 9"'' pentaquark using SRs have appeared El El 
IHlinilini- However, these calculations were either restricted to low dimension operators 
jSl El El El or they used interpolating currents which did not have the most suitable 
quantum numbers to project onto the 0"*" p, ITUj. 

The dynamics associated with the instanton, the 't Hooft interaction, has been suc- 
cessful in understanding the spectroscopy of four-quark [TTj and H-dibaryon [T21 IHl states 
and it was important for the spectroscopy of the pentaquark . Moreover, instantons are 
crucial for understanding chiral symmetry breaking in the strong interactions fT^ ITB] and 
lie at the basis of chiral soliton model for baryons, which has predicted the pentaquarks 
and their peculiar properties, e.g. small widths and masses [TTj . 

In the SR calculations thus far, the contribution from so-called direct instantons, has 
not been investigated ^. It is well known that direct instantons play an important role 
in the SRs calculations to determine the properties of the pseudoscalar mesons and the 
nucleon octet baryons [121 120] • We showed, in a mixed model-SR calculation, that they 
might be also important for the pentaquarks [21]. 

We perform a calculation for the pentaquark SRs which takes into account operators 
up to dimension d = 13 and direct instanton contributions, and leads to evidence for a 
positive parity state whose mass is close to the observed mass. 

2 The standard OPE contribution to the sum rules 

The QCD sum rule approach starts from the correlator of some relevant current, 

n(g2) = tld''x e^^-^(O|Tr/e(x)r/0(O)|O) = qU,{q'') + U^iq') . (1) 

Here ?7e represents a current with non vanishing projection onto the pentaquark state. 
We use the conventional notation q = ■ q. 

The use of the narrow resonance approximation, 

ImH(g2) = vrA|(g + Me)6{q^ - M|) + e{q' - sl)[ qlmU,{q^) + ImU^iq') ] , (2) 

where Mq is the mass of the pentaquark, Ae its residue, Sq the threshold, and the appro- 
priate dispersion relations lead to the so-called chirality even 

- ds' e~^'/*^'lmHf ^^(^2) = Xle-^l/^'' , (3) 

TT Jo 

and chirality odd 

- r° ds' e-^'/^'lmH^^^(s2) = XlMee-^'e/^'" (4) 

TT Jo 

"'^The direct instanton contribution of ref. should vanish due to the PauU principle for the quarks 
in the instanton field. 
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sum rules. 

Our choice of current in the pentaquark correlator is 

Ve = (^afgeabcebde[{u^Cde)j5Cs^][u^Cj5dg] , (5) 

whose structure corresponds to the A-state of refs. [HllHIl, which consists of the product 
of a scalar uds triquark and a pseudoscalar ud-diquaik. It can be easily seen, that this 
current has the same structure as that of ref. [H] except for the 75 in front of the strange 
quark field. The consequence of this similarity is that for the chirality odd sum rule our 
results become identical to their results, if we restrict the calculation to low dimension 
operators, take into account our different normalization and an additional negative sign 
due to negative intrinsic parity of our current, Eq. (0). 
We have also considered the current, 

vi = ^acd[{ulCjf,db + ulC-f^da)-f5l''Csl][u^C-f5dd] , (6) 

which corresponds to the B-state of refs. E] and contains a vector uds triquark 
and a scalar uci-diquark. This current has also negative intrinsic parity. However, our 
analysis has shown that the B current coupling with G"*" is weak and therefore no definitive 
conclusion about the values of the mass and residue can be drawn from the consideration 
of its correlator. We will therefore not discuss it further. 

Let us proceed with the analysis of the chirality odd SR for the pentaquark, Eq. (jH). 
This SR is directly related to the mass of the state and usually is more stable than the 
chirality even SR for the ground state baryons (201 ^^^1 triquark uds states pij . 

After Borel transforming the pentaquark correlator has dimension d = 13. The calcu- 
lation of the chirality odd sum rule will be performed taking into account operators up 
to dimension d = 13. To this order we obtain good stability in the OPE result. The op- 
erators of higher dimension {d > 13) appear in the sum rule multiplied by inverse powers 
of the Borel mass, thus, in the interesting Borel mass region, their contribution is small 
and can be safely neglected. 

With our interpolating current, the relevant trace to the chirality odd sum rule is 
expressed as 

Tr(O|T?70(x)r/e(O)|O)odd = -{if^eabcebdeeafgea'b'c'eb'd'e'ea'f'g'TT['y5CS'^;J{-x)C'y5] 
X (TT[CS2TCSi,]TT[CSfJ C755/>75] + Tt[C S:;f C St,,i^]TT[C SfJ C^,S%,] 

-TT[CS:ifCSi,CS]'JC^,S'^^n5] - Tr[CS,^rC^4,75C^^;>TC75^?e']) (7) 

where the superscripts on the quark propagator mean the quark flavor and a,b,c, are 
the color indices. In Fig. ^ the diagrams which contribute to the chirality odd SR up 
to d = 13 are shown. In order to calculate the correlator to a certain order we need 
to consider the quark propagator to the appropriate dimension. In Fig. |21 we show the 
corresponding OPE diagrams for the quark propagator which lead to 

S'abi^) = -t{0\Tq,ixMO)\0) 

= 5ab {xF^ + F^) - igGab (^^M^ + ^i^v^) 

-m^gGZo^. ( ln(-x2 AV4) + 27sm) , (8) 
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(a) (b) (c-1) {c-2) (d) 




Figure 1: Diagrams contributing to the chirality odd pentaquark sum rule in our calcu- 
lation. The quark propagator on top corresponds to the strange quark and x denotes a 
quark mass insertion. 



a, b are the color indices and g = (?c/327r^. The two functions are given by 



1 mg{qq) . niqx'^ 

H — h^ ^. ^. 9c{q(y ■ Gq), 



48 ' 27-32' 
+ ^-TTT ~ T7^9c{q(y ■ Gq) + i 



47r2a;2 ' ^ 12 192 

,2 



2 4 

29-33 



29 - 37r2 



(GV ln(-x^A74) + 27£;m - - 



where is the Euler-Mascheroni constant and we take A = 500 MeV 
for massless m, d quarks, = F/. 

Our result for the SR including operators up to dimension (i = 13 has the form 



(9) 

Note, that 
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(10) 



Each term corresponds to a diagram in Fig.E The residue is defined by Ae = (47r)^Ae. 
The contributions from the continuum are given by the following functions : 



En{M) 



T{n + l)M2'^+2 Jo 



dx e-^/^'x'^ 



(a) 






(c) 



Figure 2: (a) The terms Ff in the quark propagator, (b) the terms in the propagator, 
(c) two last terms in Eq. |S1 



with %l){n) = l + l/2 + -- - + l/(n — 1)— -"yEM- The numerical values of various quantities 
in the sum rule will be given in below. 

We have checked that our result Eq. (fTUj) taking into account only contributions of 
operators with dimensions d < 9 is in agreement with the previous calculations [HI, jH|- 
We should mention that due to specific spin and color structure of our current Eq. ^ the 
potentially important contributions from operators (iptp)"^, {iljilj)'^{iljgccr ■ Gip), ms{gc<iCF ■ 
Gq)^, and {ss){gcqcr ■ Gq)^ can not appear in chirality odd sum rules. The absence of 
the contribution of these operators can be seen by direct analysis of the terms in Eq. (7). 
We would like to mention that their appearance depends strongly on the structures of the 
interpolating current. For example, in the recent paper j2S], where another pentaquark 
current has been used, it was shown that these operators give non- vanishing contribution. 
We do not include the contributions proportional to (g^GGG) and {g^GG)"^, which in 
the OPE are related to terms of higher orders in the expansion in the strong coupling 
constant, and therefore their contributions are expected to be very small for the light 
quark systems (see for example the discussion in ref. jH] on the three-gluon condensate 
contribution). This statement is in the agreement with a general observation that pure 
gluonic operators are not very important in QCD SRs for hadrons consisting of light 
u—,d—, and s— quarks [21], |2S1- 

3 The direct instanton contribution to the sum rule 

In addition to contributions of power type, arising from the OPE expansion, there are 
exponential contributions coming from direct instantons contributions to the correla- 
tors [T^ |2ni ■ They can be calculated by using the following formula for the quark propa- 




(11) 
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gator in the instanton background in the regular gauge 

Sir\x,y) = A,{x,yh,^,{l + ^,){UT-T;^UX,, (12) 

where 

Ag{x, y) = - \Q^2^J ix - zo)<Piy - zo) 
and ^ 

(f){x - Zq) = — -r— . 

[{X - ZqY + p^\-^'^ 

Note that p stands for the instanton size and Zq the center of the instanton ; U represents 
the color orientation matrix of the instanton in S'f/(3)c and r^^y are SU{2)c matrices ; 
m* = m'^uj. — 27i'^pl{qq)/3 is the effective quark mass in the instanton vacuum and 
the current quark mass. The final result should be multiplied by a factor of two to take 
into account the anti-instanton contribution, and has to be integrated over the instanton 
density. 

To leading order in the instanton density, the direct instanton contributions arise from 
two body ud, us, ds and three body uds quark zero mode propagators in the correlator 
Eq.(IT)), as shown in Fig. 3. 




Figure 3: An example of instanton two- and three-body contributions to the correlator of 
the pentaquark currents. 

The final result for two body instanton contribution is 

MM) = -^?f4%^5[/6(g)], (13) 

where Shuryak's instanton liquid model for QCD vacuum with density n(p) = neffS{p—pc) 
has been used and B[fQ{Q)] is the Borel transform of fei^Q) which is defined by 

J ^'" x^[.g+pg]3[(x-.o)^+pg]3 ' 

where pc is the average instanton size. There are two types of singularities in Eq. p4j) . 
One of them is related to the pole at the origin = 0, the other is due to the pole 
at finite distance from origin ~ —p^- The pole at = produces, after Fourier 
transforming, power terms in in addition to the exponential type direct instanton 

contributions exp(— Qpc), arising from finite distances. One should carefully subtract 
that contribution to avoid double counting with the standard OPE terms. We follow the 



5 



procedure suggested in ref. |2d] for the analysis of direct instanton contributions to heavy 
quark decay. More specifically, for a general integral 

where S{x) contains no singularities for complex x^, we use Feynman's parameterization 

T{a)T{j3) J ' Jo (t(l - t)x2 + zl + p2)-+/3 ' ^'''> 

where Zi = Zq — tx, but consider only the contribution from the pole at 

= ~{zl + pIMI - t)]-\ (17) 
The Borel transform of the function /g is given by 



+ 10X^ + 5X'' + X'j , (18) 

where X = {1 — t)/t. Note that only the contribution from the pole at finite quark 
separation has been considered. 

We have also performed the calculation of the three body contributions induced by 
instantons, Fig. 3b, and found that they vanish as the result of the cancellation between 
diagrams with different uds combinations in the pentaquark current. The two body ud 
instanton contribution also cancels. The only non vanishing contributions arise from the 
two body us and ds terms of Eq. p3|) . The behaviour of the different contributions is 
associated to the Dirac structure of our current Eq. which includes both scalar and 
pseudoscalar ud diquarks with equal weights. The instanton induced contribution is very 
sensitive to the parity of the state [2Ij and it flips sign when the parity of the state changes. 
Of course, the three-body instanton induced forces might give non-zero contribution for 
other choices of pentaquark currents, for example, for currents with derivatives. 

We should mention that three body instanton terms induce forces which give non 
zero contributions to the mass of some specific triquark uds pentaquark clusters ^1 1^ 
and furthermore, they give non vanishing contribution to the G mass within the bag 
model 



4 Numerical analysis 

We use the following values for parameters at the normalization point 2 GeV [8| (see also 
recent discussion about uncertainties in values of various condensates in j^Hl) 



{uu) = -(243 MeV)3 ^ 

b = {glG^) = 0.88 GeV^ 

igc{ua ■ Gu) = ml{uu) = O.SGeV'^ {uu) , 

{ss) {sa ■ Gs) 

{uu) {ua ■ Gu) 
= 111 MeV, 



/. = 0.8, 



(19) 
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Figure 4: The mass of the pentaquark 
obtained without direct instanton con- 
tributions as a function of the Borel pa- 
rameter in different orders of the OPE 
expansion. 
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Figure 5: The residue of the pen- 
taquark obtained without direct instan- 
ton contribution as a function of the 
Borel parameter in different orders of 
OPE expansion. 



with pc = l.Q GeV ^ for the average instanton size in the QCD vacuum. In the numerical 
estimate of the direct instanton contribution the relations of the instanton liquid model 

ISO! 

nrr ?. 777* 1 

(20) 



Ueff _ 3 rUg _ 1 



are used. 

In Figs. 4 and 5 the result of the calculation of the pentaquark mass and residue within 
the standard OPE expansion for the different orders in operator dimensions is shown. In 
Figs. 6 and 7 the mass and residue of the pentaquark as a function of the value of the 
Borel parameter with direct instanton contributions are shown. In Fig. 8 we present the 
results of the calculation of the OPE and the direct instanton contributions to the left- 
hand side of the SR, Eq.(j3]). All curves are given for a value of the threshold sq = 2 
GeV. We chose this value of threshold because the stability was best. From the fit of the 
sum rules we arrive at the following values for pentaquark mass: Mq+ = 1.66 GeV for 
d = 7, Me+ = 1.75GeV for d = 9, Me+ = 1.73GeV for d = 11, and Me+ = 1.75 GeV for 
ci= 13 2. 

One important result of our calculation is in the change of the sign of the squared of 
the residue when increasing the dimension of the operators which contribute to the OPE. 
Thus, for d = 5 the sign is positive, while it becomes negative for higher dimensions. In 
particular, the contribution from the dimension d = 7 operators is crucial for inverting the 
sign. Due to the negative intrinsic parity of our current Eq. (0), the negative (positive) sign 
of the squared of the residue implies positive (negative) parity for the state. Therefore, 

^Once the instanton contribution is included the stabiHty in the Borel parameter for the SR up to 
d = 5 operators disappears (see Fig. 6). Therefore it is not possible to extract the value of pentaquark 
mass from the SR with only up to d = 5 operators. 
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Figure 6: The mass of the pentaquark ob- 
tained with direct instanton contributions as 
a function of the Borel parameter for differ- 
ent orders of OPE expansion. 
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Figure 7: The residue of the pentaquark ob- 
tained with direct instanton contribution as 
a function of the Borel parameter in different 
orders of the OPE expansion. 



our final result for the residue presented in Fig. 7 shows that one can arrive to the wrong 
conclusion about the parity of the pentaquark state [Hj, if one takes the decision based on 
only the contributions of low dimension {d = 5) operators. We also stress the necessity to 
include high dimension operators to get good convergence for the sum rule. It is evident 
that this effect is directly related to the high dimension of the pentaquark current. 

Once we include the contributions from the high dimension operators and the instan- 
tons, our result for the B"*" pentaquark mass, Mq+ ^1.75 GeV, is higher than was given 
by previous SR calculations jHl El El but still in rough agreement with the available ex- 
perimental data, if one admits about 10% accuracy in the predictions of the SR approach 
due to uncertainties in the values of the various condensates, the mass of the strange 
quark, the contribution from higher dimension operators d > 13, higher order pQCD 
corrections, etc. Furthermore, some additional effects such as the mixing between various 
pentaquark states [1311211, which are beyond the scope of the present paper, might give 
some additional contribution to the mass of the G"*". 

We also note that in our calculation the pentaquark has positive parity in agreement 
with the soliton model prediction ^71. Our estimate for direct instanton contribution 
is done within Shuryak's instanton liquid model. We have found that the instanton 
contribution for the full SR is rather small, but can give a large contribution to it when 
one considers operators only up to dimension d = 5 (see Figs. 4,6). The smallness of 
the instanton contribution to the full SR is mainly related to the large mass of the Borel 
parameter M ^ 1.7 GeV, where we obtain the plateau of stability (Fig. 6). In this region 
the instanton contribution is small in comparison with the contribution from the high 
dimension operators in the OPE (Fig. 8). 

There is a significant dependence of our results on the value of threshold. This is 
a common feature in all the studies about the properties of the pentaquark within the 
QCD sum rule approach. In our case, we have chosen Sq = 2 GeV to satisfy the physical 
requirement of having a large stability plateau. 

In summary, we have shown the analysis of the QCD sum rules for the pentaquark 
current including high dimension operators in the OPE and direct instanton contribu- 
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Figure 8: The different OPE orders and direct instanton contributions to the left-hand 
side of chirahty odd SR . 

tions. Our results conclude that the role of the high dimension operators is important 
for obtaining a positive parity for pentaquark state. Our calculation though produces 
a bound state whose mass is higher than the experimental observation. More sophisti- 
cated models and probably states mixing [HI might reduce the obtained value to the 
observed one. 
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